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Abstract 


This article is a continuation and an extension to the previous work covered 
in Ref. (1) that deals with the thermal behavior of a tubular solar collector. 
A special analytical solution is provided for the time-wise response of the 
circulating fluid temperatures when a sudden step change of the input solar 
radiation is Imposed and remains constant thereafter. An example which demon 
strates the transient temperatures at the exit section of a single collector 
with two different flow patterns is presented. This study is used to supple- 
ment some numerical solutions to provide a fairly complete coverage for 
this type of solar collector. 
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The Transient Thermal Response of a Tubular Solar Collector 


Introduction 

The study of transition states in the performance of solar collectors is 
becoming increasingly important with the new introduction of solar module con- 
trol techniques. The temperature control versus the circulating fluid rate, 
for example, is directly related to the transient behavior of the module at 
specific input conditions of solar irradiancy, ambient and inlet fluid temp- 
eratures. 

In order to study the transient performance of any heat exchanger, a 
sufficiently general analytical model must be established which is both an ade- 
quate idealization of the physical system and capablr of reasonably simple 
mathematical description. There are basically three common methods in the 
literature for solving the transient performance problem. These are classified 
as : an analog computer method, a numerical method, and an analytical method. 

The choice of the method to be used depends on the complexity and size of the 
problem, the accuracy required, the means available, and the experience of the 
investigator. First, analog computer methods (Ref. 3, 4 & 5) provide 
a realistic and accurate analog model of the physical system, but require 
a significant amount of advanced electrical circuitry. Second, numerical 
solutions employing the practical method of finite differences has become in- 
creasingly popular with the availability of high speed digital computers and 
its associated technology. This technique is very useful for almost any heat 
exchanger problem, and is the only method that can be used if a high degree of 
accuracy is required or if the problem at hand is of considerable size and com- 
plexity. One using this method has to solve algebraic equations (nodal equations) 
instead of the partial differential equation. The nodal equations can be derived 
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by purely mathematical methods or by energy considerations. Third, analytical 
methods, as evident from Ref. 2, present in general a substantial mathematical 
complexity. However, a special case has been found for tubular solar collectors 
which can be handled by standard orocedures and this is the subject of this 
article. The analytic solution of the transient behavior can later be used 
to supplement the numerical solutions so as to provide adequate coverage for 
this type of solar collector. 

Collector Description and Flow Patterns 

The tubular collector, as shown schematically in Figures 1 and 2, is com- 
posed of three concentric tubes*, an inner tube, an absorber tube, and a cover 
tube. The annulus space between the absorber and cover tubes is highly evacu- 
ated to minimiza convection and conduction losses. The absorber tube surface is 
coated with a selective material to reduce the outward long-wave radiation losses. 
In flow pattern (1), as shown in Figure 1, the circulating fluid starts from the 
inlet section of the inner tube. At the closed end of the collector, the fluid 
reverses its direction and passes in the annulus spacing between the inner and 
absorber tubes. In flow pattern (2), as shown in Figure 2, the fluid path is 
reversed from the abov?\ Both flow patterns alternate in each collector module 
as shown in Figure 3. For more irradiancy augmentation, the set of collectors is 
mounted with lateral spaces separating them from each other with a highly reflec- 
tive back reflector. 

Assumptions and Analysis 

It is widely recognized that the analysis of the transition states of heat 
exchangers does not respond to an easy mathematical treatment, and that simpli- 
fying assumptions are necessary in order to achieve any type of solution. 
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The following idealizations are made in the model to yield a transient 
solution with an accuracy adequate for engineering purposes: 

1 - The heat flux and the temperature of each of the fluids and the walls 

are functions of both time (9) and axial distance (X) from the collec- 
tor inlet. This makes the problem one dimensional. 

2 - The densities and specific heats are constants throughout the system. 

3 - The thermal conduction of the tubes' material is assumed zero in the 

axial direction parallel to the flow and is considered infinite in the 
radial direction. 

4 - The heat transfer coefficient between any tube surface and its adjacent 

fluid is uniform and constant over the axial direction for a fixed 
fluid flow rate. 

5 - The heat transferred at any section due to the fluid thermal conduc- 

tivity in the axial direction, within the fluid, is negligible compared 
to other heat fluxes from or to the fluid at the same section. 

6 - The heat capacitance of the three tubes of the collector, namely the 

inner, the absorber and the cover tube, are negligible relative to the 
heat capacity of the fluid. This is a good simplifying assumption if 
a liquid is used as the circulating fluid in the collector. In case of 
using gases or air as the solar collector medium, this assumption is not 
recommended since the wall capacitance has large effects on the thermal 
lag. 

7 - Initially, steady state conditions were prevailing before a sudden input 

variation in the solar radiation was superimposed. 

8 - The fluici flow rate is constant throughout the collector. 

9 _ passes (1) and (2) of the inner and outer fluid, respectively, are of 

equal flow area; thus equal flow velocities in each pass are assumed. 
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Virtually all of these idealizations are also required for the 
analysis of the steady state performance. 

Based upon these idealizations, the differential equations relating the 
system temperatures are derived from energy balance and heat transfer rate 
equations applied to a segment of the collector of length dx. The derivation 
of these' differential equations is presented in Appendix A. The incident 
solar radiation is then applied in the form of an input step function. The 
transient problem for both flow patterns is then mapped on the frequency domain 
instead of the. time domain by using Laplace transformation as explained in 
Appendix B. In Appendix C, the inverse Laplace transformation is presented using 
the Residue theorem in complex number algebra, for one of the two flow patterns, 
namely pattern (1). A similar approach is followed for flow pattern (2) and is 
presented separately in Appendix D. 

Numerical Example 

The following numerical example will show the difference in the thermal 
response between the two flow patterns of the tubular solar collector . The 
data were arbitrarily abstracted to be as close as possible to actual running 
conditions. However, the conclusions may be generalized at any other relevant 
conditions. 

Starting with inlet water temperature 
Ti (0,0) = 70°c for flow pattern (1) 

or 

To(0,0) = 70°c for flow pattern (2) 
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Ambient temperature , 


* * 


30°c 


Tube length 1.067m 

Absorber tube diameter . . 0.041m 
Cover tube transmissivity 0.91 
Absorber tube absorptivity 0.85 

Mass flow rate 5 Kg/hr 

Water specific heat .... 11.63 x lO^Kwhr/Kg°K • 

Sky temperature 4°c 

l eeder tube diameter . . . 0.029m 
Cover tube diameter .... 0.051m 
Cover tube emlsslvlty ... 0.90 
Absorber tube emlsslvlty. . 0.1 

Wind speed • 11 Km/hr 

Augmentation factor .... 1.64 

Irradiation intensity varies in the form of a step function and is given by 
0 .0 < 0 
0.75 Kw/m 2 0 <0 < w 
The relevant constants K-j , Kg, K 4 (0) were calculated from Eqs, A-20 & A-23 using 
the heat transfer coefficients data (Ref. 6); 

0 < 0 °K/m 
O£0<« °K/m 

AK 4 (0=0) = 8.3537 °K/m 

C = (VM = 0.0085495 
— 2 

Vi=V 0 =V=7. 5698 m/hr 

Since the initial conditions were at steady state. Equation B-9 was solved as 
in Ref. (1) and gave the difference ['Tq( 0,0) - T^ (0,0) ] as -0.818 °K 


K 1 = 0.8853 m" 1 K 4 =[5.0869 

Kg = 0.9024 m" 1 I 13.4406 
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The first six vectors Z in Equation C-17 and the corresponding poles p are 


calculated as: 

Z = i 1 .8274 j 
Z = 2.1065 + 7.5750 j , 
Z = 2.7063 + 13.9438 j , 
Z = 3.0754 +20.2691 j , 
Z = 3.3437 + 26.5781 j, , 
Z = 3.5548 + 32.8788 j , 


p = -5 ,0>55 

p = -14.4634 + 26.0865 j 
p = -16.4525 T 49.0227 j 
p = -17.7222 + 71.5935 j 
p = -18.6569 + 94.0439 j 
p = -19.3968 + 116.4385 j 


The first six poles constitute 97% of the infinite series sum required for flow 
pattern (1) and are sufficient for temperature calculations. For flow pattern (2) 
many more transition terms are required for calculations in the vicinity of the 
(9 = 0) region. These extra terms will have a very rapid diminishing effect as 
0^ 0.1 hr . 

The temperature T o (0,9), from Equation C-16, for flow pattern (1) is given 


by: 

T o(°,9) = T Q (0,0) + 8.7826 - 6.5290 


-5.0655 9 
e 


-14.4634 0 

+2.7288 e . cos(2.1593 + 26.0865 0) 


-16.4525 9 

+0.4480 e ’ . cos(l. 5462 - 49.0227 0) 


-17.7222 9 

+0.5363 e ’ . cos (2. 3276 + 71.5935 9) 


-18.6569 0 

+0.1964 e cos(l .2973 - 94.0439 9) 


-19.3968 9 

+0.2447 e cos (2. 3688 + 116.4385 0) 
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and that for flow pattern (2), from Equation D- 9,. is given by: 

-5.0655 0 

T i (0,0) = T.(0,0) + 8.7826 - 14.6569 e 


-14.4634 0 

+2.9974 e * cos(l .7492 - 26.0865 9) 


-16.4525 0 

+1 .0818 e cos (2. 0444 - 49.0227 0) 


-17.7222 0 

+1.0681 e res (1.6352 - 71.5935 0) 


-18.6569 9 

+0.6108 e cos (1 .8797 - 94.0439 0) 


-19.3968 0 

+0.6429 e cos. (1.6005 - 116.4385 0) 

The outlet temperatures are plotted as shown in Fig. 4 for both flow patterns. The 
trend of the two outlet temperature curves is the same for any other conditions of 
flow rate, ambient temperature or solar irradiancy. 
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Summary and Conclusions 

To summarize the major features of this work the following conclusions 
are made: 

1 - A pure analytic procedure using Laplace transform has been established 

for the transient thermal solution of the tubular solar collector with 
two different flow patterns. The input solar flux is introduced in the 
form of a step function and the resulting collector response is pre- 
sented in each case. The numerical results of a selected example are 
plotted in Fig. 4 for comparison. 

2 - Comparison of the exit temperature from the collector. Figure 4, indicates 

the distinct difference in responding to a sudden change of solar radiation. 
The outlet fluid temperature of flow pattern (1) responds by a much faster 
rate in the early stages since the fluid leaves the collector as soon as 
it absorbs the useful energy from the annulus area facing the sun. On 
the other hand, the fluid temperature of flow pattern (2) lags behind by 
the "residence" time. The latter is the time elapsed for an element of 
fluid to pass through the collector from the inlet section to the end 
section. The residence time 9 $ is given by 



and is equal to 8.46 minutes for this example. This means that approx- 
imately 8 minutes have to elapse to discharge the relatively colder 
fluid in the center tube before a significant temperature gain is ob- 
served for flow pattern (2). 

3 - The steady state collector analysis, described in Ref. (1), leads to 
the conclusion that both flow patterns, though having different tempera- 
ture distribution inside each tube, will eventually end up with the same 
steady state outlet fluid temperature. Since the temperature growth for 
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flow pattern (2) is lagging behind that of flow pattern (1) and catch- 
ing up at some infinite time, the rate of temperature increase for flow 
pattern (2), at the end of the "residence" time has to be much faster 
than the early rate of response for flow pattern (1). 

4 - The steady state time registered experimentally by the collector manu- 

facturer and by the solar simulator team at NASA Lewis Research Center 
for this collector type, ranged from 30 to 60 minutes according to flow 
rates, ambient and solar flux conditions. This relatively long steady 
state time is due to the fact that 65% of the collector volume Is filled 
with liquid fluid that hinders a fast response. The analytic results 
of the given example seem to agree with the above finding; the outlet 
fluid temperature reached 97% and 92.5% of its steady state value for 
flow patterns (1) and (2), respectively after 40 minutes of exposure to 
sudden solar radiation. 

5 - The series solution presented in this work depends entirely on the para- 

meters K-| , C , L and V which in turn depends on the physical dimensions of 
the collector, the heat transfer coefficients and the fluid flow rate. 

With fixed collector dimensions, the higher the rate -of fluid flowing in, 
the higher the conductance coefficients and the faster the response to 
solar flux changes will be^even though the useful temperature rise across 
the collector will decrease. On the other hand, with a fixed fluid flow 
rate, the smaller the area of the flow passage and the collector size, the 
higher the fluid velocity, the higher the conductance coefficients and 
the faster the rate of heating will be^ since a smaller amount of liquid 
is residing inside the collector. 
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6 - The analytic solution can now be applied to different Input conditions 
of flow rate, solar flux, ambient temperature or collector dimensions 
for optimizing the collector performance. Moreover, other step changes 
of the solar flux, such as sudden shading due to clouds passing by and 
their effects on the collector cooling rate, could be followed by the 
same expressions derived herein. 

In conclusion, it should be pointed out that the analytic solutions obtained 
in this article could be handled with a digital computer to give adequate cov- 
erage for this type of solar collector. 


10 
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Nomenclatures 


cs 


C f 

D 

F 

I 

4 

L 

m 

P 

r 

r 

T 

T 


sky 

'"amb. 


ac 


ao 


'cv 


U . 
01 


rd 


Kwhr/Kg°c 

Kwhr/Kg°c 

m 

> 1 
Kw/m 2 

m 

Kg/ hr 


cross sectional area m 

specific heat of tube material 
specific heat of flowing fluid. 

Tube diameter 

Augmented Radiation Factor 
Direct incident solar radiation 

tube length 
fluid mass flow rate 
Laplacian frequency 
reflectivity 

Transformed temperature on frequency domain 

Temperature °K 

Sky temperature °K 

Ambient temperature °K 

Radiative heat transfer coefficient between the absorber and the 

p 

cover tubes Kw/m °c 

Overall heat transfer coefficient between the absorber tube and 
the outer fluid in the annulus 
convective heat transfer coefficient between the cover tube and 
the ambient air 

overall heat transfer coefficient between the outer fluid and 
the tnner fluid Kw/m 2 °c 

Radiative heat transfer coefficient between the cover tube 
and sky. 

Fluid velocity. m/hr 

Distance along the collector tube m 


Kw/m 2 °c 


Kw/m 2 °c 


Kw/m Zo c 
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a 

absorptivity 

£ 

emissivity 

r 

transmissivity of cover tube 

S 

tube thickness 

cr 

Stefan Boltzmann constant 

P 

density 

Subscripts 

a 

absorber 

e 

cover tube 

i 

inner fluid, passCl ) 

0 

outer fluid, pass (2) 
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AMBIENT TEMPERATURE 
SKY TEMPERATURE 



Fig, la. Collector configuration 
(flow pattern 1) 



o distance along TUBE * L 


Fig. lb. Temperature distribution 
along one collector (flow pattern 1) 
at time 0 
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ambient temperature 

SKY TEMPERATURE 



Fig. 2a. Collector configuration (flow 
pattern 2) 



Fig. 2b. Temperature distribution 
along one collector (flow pattern 2) 
at time 0 
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Fig. 4. Transient response of the outlet 
fluid temperature of the tubular solar 
collector as derived from the analysis 


IS 
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Appendix A 


Derivation of the system differential Equations for the Transient Solution 

For a segment of the collector of length dx at the time Interval from 9 to 
{ 0+dG ) , the rate of heat flux is divided, as shown in Figure A-l» whereby 
dQ-| * The total radiation on the cover tube from all sides including 

Irradiancies from the back reflector and reflections from adjacent 
collectors: 

dQ 1 = I , F , (D c dx) . . . (A-l ) 

dQ 2 = The total outward reflection loss from all sides of the cover 
tube, 

dQ 3 = the energy absorbed by the cover tube including absorption from 
multiple reflections within the evacuated space 

d( b = a c + W r * d ^l * • • ^ - 2) 

c a 

dQ^ = the energy absorbed by the absorber tube including absorption from 

multiple reflections within the evacuated space, 

% ** r" £ > * d Ql . • • • CA-3) 

c a 

dQ 5 = the long-wave radiation exchange between the cover and the absorber 
tubes. This is absorbed by the cover tube almost entirely. 

. . (A-4) 


. . . (A-5) 


. . . (A-6) 


where 


d <)5 ■ U ac • < 7rD a' dX > ^ a - T c> 


U__ = cr . f 
ac 


j4 t 4 

( V T c ) 

ac '(V T c) • 


i = i + -1 (i . ! ) 

e ac «. Vt ’ 


dQg = the sensible heat gain by the absorber tube material 
dQ fi = tt D . 8 . dx . P . C(iZa) 

0 a a a a 0 q 
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dQ 7 = The heat transfer to the outer fluid In the annulus pass from 


the absorber tube surface 


= U aQ . (7rD a dx) . (T a -T 0 ) . . . (A-7) 

dQ 8 * The combined radiation and convection loss from the cover tube 
to the ambient air and sky. 

dq 8 = (7rD c -dx) . [u cv (T c .T amb ) ♦ U rd (T c -T 5ky )] _ . . (A „ 8 ) 


where 


t 4 t 4 

u rd ■ £r ' « e -U->£ 

' c " ‘sky 


dQg = The sensible heat gain by the cover tube material. 

dQ 0 = S„.dx. P n . C_ / 9T C\ 

^9 c c c c ) 


. . . { A-8 ) 


. . . (A-9) 


. . . (A-10) 


dQ^ = The sensible heat gain by the outer fluid In the annulus pass. 
This is determined by the "Lagrangian" approach following the 
elementary fluid mass as it passes by the section bounded by X 
and (X + dx) during the time Interval from 9 to {0 + d9) 


dq 10 s [ A cs, 2 P f ,C f (*|^-) + * C f (^-)] 


. . .(A-ll) 


where the (-) sign is for flow pattern (1) and the (+) sign for flow 
pattern (2). 

dQ^-| = The heat transfer to the inner fluid, 

dQ-j-j = U 2i (tto^xM^-T.) . . . (A-12) 

dQi g = The sensible heat gain by the inner fluid, pass (1), This is 

determined, as in dQ^ Q , by following the elementary fluid mass as 
it passes by the section bounded by X and (X + dx) during the time 
interval from 0 to (9 + d9), 


= [ A cs.r p f ,c r hfe 1 ^ - 1,1 


18 
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where the (+) sign Is for flow pattern (1) and the (-) sign for flow 

pattern (2). 

and 

dQ^o * The sensible heat gain by the Inner tube material. 

Applying the first law of thermodynamics to the system componets will 
yield: 

dQ 3 + dQg - dQg - dQg s 0 . . .(A-14) 

for the cover tube, 

d(fy ~ dQg - dQg - dQy * 0 , . . (A-15) 

for the absorber tube, 

dQy - dQ-jQ - dQ-| i s 0 . . .(A-16) 

for the outer fluid, and 

dQ^ - dQi g - dQ-| 3 * 0 . . .(A-17) 

for the Inner fluid. 

In this work, it has been assumed, without great loss of accuracy, that the 
heat capacitance of the three tubes of the collector are negligible relative to 
that of the liquid fluid. Accordingly, the number of differential equations 
decreases to two instead of four since 

dQg = dQg = dQ.| 3 = zero . . . (A-18) 

and equations A-14 and A-15 are reduced to: 

CC a c + [ - yjr - ) .I.F.(P C dx)] + [(7TD a dx).U ac .(T a - T c ) ] 

51 C 

-7TD c dx [O c -CT c - T am(j ) + U rd CT c - 0 

' I - F • C° c dx ) ] - [(7rD a dx) . U ac (T a -T c ) ] - 

- [(’r0 a dx).U ao .( V T o )] - 0 


OPL Technical Memorandum 33-781 


19 



These can be written as: 


T c * » + b T a 

T c ” fT, - (f-1 )T 0 - g 
where b* f & g are characteristic constants given by: 
a c r a Tc 

L'~c T 
a = 


[<“c + + ( U cv .T amb + U rd T sky) j 

f ii _i_ n a . /D~ 1 1 \ 1 


[ U cv + U rd + U • u ac>] 

C 


b = 


(?• ' U a=) 

c 

< u cv + u rd + ^ • u ac > 
u c 


• • » (A-1 9 ) 


. . . (A-20) 


u an 

f = (1 + Tj ) , 

ac 


g . LI 


a a T c 


w(U a „. u a ) \ l-r,r„ 
' ac fp ' \ a c 

c 


On the other hand, Equations A-16 and A-l 7 will become the twr governing 
differential equations written as 

[U ao . ( W D a dx ) . CT a -T 0 ) ] - [U of ( 7T 0i dx) (T,,-^ ) ] - 


" I a cs .2 p f c fW * * c f a x i dx “ 0 
«. 

where the (-) sign is for flow pattern (1) and the (+) sign for flow pattern (2) ; 

[ U oi (7TDjdx)(T 0 -Tj ) ] - [A csl fi j, C fil“3g^) ± * C f^"gX’^ dx = 0 

where the (+) sign is for flow pattern (1) and the (+) sign for flow pattern (2). 

The continuity (mass conservation) equation Is expressed as 

^ = P f A csl ,V i “ P f A cs2 ,V 0 ' ’ •[&”%*) 
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;■ id together with Equations A-19 and A-20, the simultaneous partial differential 
equations governing the thermal behavior of the tubular solar collector ere 
expressed as: 


3Tj 1 3T/ 

inr + T, mf + K 1 T 1 - ¥o ■ 0 

aT. 1 aT. 

±sr - V 0 TT + K 1 T 1 • k 3 t o 'V 0 


where 


K 3 = K 1 + h 
K 4 = K 2 


h = 


ttD U 
a ao 

rti C- 


. . . (A-22) 


. . . (A-23) 


the (+) sign is for flow pattern (1) and the (-) sign is for flow pattern (2). 
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Fig. A-l. Sankey Diagram for the 
Solar Collector 
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Appendix B 

Transformation of the differential equations from the time domain to the 
frequency domain. 

The partial differential equations that govern the collector inner and 
outer fluid temperatures T^ (X, 0) and T Q (X, 6), respectively, are expressed as 


+ 


3T, , 3T, 

JT V Tff 

3 T o 1 ® T o 

~ST ‘ V IT 


+ K 1 T 1 - K 1 T o ’ 0 


+ K, Ti - Kj T 0 + K 4 - 0 


. . .(B-l) 


where K-j , Kg, & are characteristic parameters, and 

V i = V Q = V . . . (B-2) 

according to the ninth assumption in the main text. The (+) sign is for flow 
pattern (1) and the (-) sign for flow pattern (2). At steady state conditions, 
the above set of equations are reduced to the same set derived in Ref. (1), 
namely 


+ Ki T 4 - K, T 0 = 0 


'1 


. d 1 
i 33T 

.T ft 

+ + K n T„. - K 0 T^ + K/, = 0 


. . . (B-3) 


'1 'i ,x 3 'o ' 4 

where the (+) sign is for flow pattern (1) and the (-) sign for flow pattern (2). 

The following Laplace transforms will convert the set of equations in (B-l) 
from the time 0 domain to the frequency p domain to ease their algebraic manipu- 
lation. These are 

00 


^ (X, P) = / 


00 


T_ (X, P) = / 


-p0 

e T i (X, 0) d0 


-p0 

e T q (X, 0) d0 


. . (B-4) 
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Multiplying equation (fi-1) by e “ p0 and integrating from (0=0) to (0=«) 
gives: 


+ 9 ^i / p - *” *1" » ( X , 0) 

1 sr + (K i + f> T i - K i T o - V- ■ 


. . , (B-5) 


+ 3 T o 


“ K d TJX, 0) 


-;-r- (K 3 J )T o + K i T i + r + 


= 0 


where the (+) sign is for flow pattern (1) and the (-) sign for flow pattern (2). 
The two simultaneous partial differential equations could be separated for either 
variable by using Cramer's rule. These are expressed as 


3 

ax 


J 7 lK3 ' K,) |r - KjXKaJ) -K, 2 ] T, + K,(piJ 


K. T n (X,0) 


v— > 


and 


rTj (X,0) p 1 _ , 8 T. 

+ [ . (Kg + y^j + y* (X,0) = 0 • • • (B“6) 

3 * (Krl<l) |r- [< K i + fx K 3 + f>- K i 2 ]To 


_a 
a x 


= -K, . 

. . . (B-7) 

where the (-) sign is for flow pattern (1) and the (+) sign for flow pattern (2). 
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Distribution Along the Tube at 9 = 0 
This fctep has to be determined first before a general solution can be made. 
Generally speaking, the Initial conditions are arbitrarily preassigned based on 
past experience. However, in this work there is no need for any initial estimates 
or guesses since we have already stated in the seventh assumption in the main 
text that Initial conditions are at steady state. Accordingly, using equation 
(B-3) , the initial temperatures ( X, 0) and T Q (X,0) are related by: 

8L(X,0) 

t / v n \ i/ i* / v n \ ~ n 


aUx.o) 


+ K 1 T 1 (X,0) - K 1 T 0 (X,0) = 0 


+ K ] T. (X,0) - K 3 T q (X,0) + K 4 (0) = 0 


. . . (B-8) 


where the (+) sign is for flow pattern (1) and the (-) sign for flow patterns (2). 
Equation B-8 can be further separated for one of the temperatures T^(X,0) or 
T 0 (X,0) as: 

d 2 T- (or T.) dT,(or T_) 

— ^ — - + (K 3 - K-| ) 2. . K] (Kg -K-j ) - K 4 {0). . . (B-9) 

where the (-) sign is for flow pattern (1) and the {+) sign for flow pattern (2). 

For each flow pattern, equation B-9 is solved subject to the boundary conditions 

at the end (X=L) as reported in Ref. (1). The result is an explicit relationship 

for the quantities T_. (X,0), T„(X,0), dT i( x »°) and dT 0 ^’°^ with the distance X 

1 0 “lx "lx 

along the tube. 
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11 - General Solution in the p-plane : 

Our Interest in this work is to study the thermal response when a sudden 
step change in the Insolation "I" is imposed,, Figures B-l and B-2 illustrate 
the resulting change in the parameter K 4 (0). The latter can be written as 

K 4 {0) = K 4 (0) + AK 4 . u (0) . . . (B-10) 

where u(0) is the unit step function at 0 = 0. 

Substituting the initial conditions given by equation (B48) in equations 
(B-6) and (B-7), will give 


ST i T 


7 + < K 3 ' MaTT ~[ (K 1 + T^ K 3 + £> ' K l 2 ] T 1 
M; t, (X,0) p> 

' - V- - V- < K 1 + K 3 + f> 


and 


. . .(B-ll) 


i. r 

d T o v 


5 t 


7 + < k 3 - K, - [tit, + f)CK 3 + f) " T c 


- T » (X ’ 0) (I- I r I "1 - 

-p— - — v — lK l + k 3 + V " T 


P\ AK 4 


. . . (8-12) 


where the (-) sign is for flow pattern (1) and the (+) sign for flow pattern (2), 
The particular "integral” or solution can thus be determined from equations 
(B-ll) and (B-12) as 


T.(X,0) _ K 1 AK 4 V 2 
~ P ” + p - (“p+r-j J rp+r^) j 

for T^ (X,p) and both flow patterns, 

and for T Q (X,p) the particular "integral” is given by 

T o (X,0) _ ^ AK 4 V 2 AK 4 V 

— P + p(p+r-j)(p+r 2 ) + (p+r^lp+rg) 
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for both How patterns where 


[< K i + + ^ “ k i 2 ] = (P +r ])(P +r 2^ 

is written for simplifying the inverse transformation 
and 


. (B-13) 


r } r z = 


(r-|+r 2 ) = 2 (C+K-j )V 


u ~y~ 


r 1 = cv + + rye 2 + K-j 2 


r 2 = CV + K-j V - V\/C 2 + 1C, 2 


• (B-l 4 ) 


On the other hand, the roots of the "characteristic" equation of the image 
differential equations (B-l 1 } and (B-12) are such that 


s l,2 + S 1 ,2 “ [( k i + 5 , H k 3 + 7’) “ K 1 ] " 0 


. . . { B-l 5 ) 


where the (-) sign is for flow pattern (1) and the (+) sign for flow pattern (2). 
The roots S-j(p) and S 2 (p) are further written in the form: 

S-j (p) = + {C + R) 

$ 2 (p) = + (C - R) 


R = 


. . . ( B-l 6 ) 


V(v +C )(7 + ct *i) 

where the (+) sign is for flow pattern (1) and the (-) sign for flow pattern (2). 
The general solution for T^(X,p) and T Q (X,p) is given by: 

S,X S-X T i (x,o) k,v 2 ak. 

TiOC.pl = A^) + A 2 tp) e + -u— + ptp +ri) Tp - r g ) ...(B-17) 
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n 


5.X FpX T (X,0) K,V 2 AK d AK 4 V 

T 0 (X,p) - A 3 (p) e + A 4 (p) e + — + pjp+r^ Up+r 2 ') + Tp+^V^) 

. . . (B-18) 

where A-j (p ) , A 2 (p) , A 3 (p) and A 4 (p) are arbitrary functions of the frequency p. 
Substituting in any one of the image differential equations such as equation 
B-5, and comparing the coefficients of exponential and absolute terms: 


A-i (p) p 

A 3 (p) « («! + C + $•+ R) 

Ap(p) p 

A 4 (p) = (K, + C + £ - R) . J 

for both flow patterns. 

Combining equations B-16, B-17, B-18 and B-19 gives 


.(B-19) 


5.X S ? X T.(X,0) KtAK/ 

T i (X,p) = A-| (p) e +A 2 (p)e^ + -^ + 


piP+r-j Hp + r 2 ) 


. . .(B-20) 


. A-i(p) p S 1 X A 2 (p) p 

Tq(Xi^J ~ (Kj+C + y+ R)® + £ (K-| + C + y 

T o (X,0) K -j AK 4 V 2 AK 4 V 

+ “p + p^p+r 1 7lp+r 2 ; + (p+^Mp+iy 

for both flow patterns except for the signs of S*-j and S^. 


S 9 X 
R) e 2 


. . . (B-21 ) 
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1(0) 


JL 


0 TIME, B 

Fig. B-l. Step change in the solar 
radiation at 0 = 0 



Fig. B-2. Step change in the parameter 
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Appendix C 

Inverse Transformation from the Frequency Domain to the Time Domain for 
Flow Pattern (1) . 

i - Boundary and Initial Conditions : 

With the image temperatures Tj(X»p) and T Q (X,p) as determined from 
Appendix B, equations B-20 and B-21 are subject to the following two conditions: 
1 , A boundary condition : 

At X = L and at all times: 

T i (L»9) = T o (L,0) 


T, (L.p) = r 0 (L,p) . . .(C-l) 

2. An initial condition : 

At the inlet section (X=0) and at all times*. the temperature T^(O,0) is 
kept fixed, i .e. , 

(0,9) = T.j(0,0) = constant 


or 


, T. (0,0) 

Ti(O.p) « 


• . (C-2) 


The arbitrary constants A-j(p) and A 2 (p) for flow pattern (1) are then giVen by: 


K 1 AK 4 .V 


[( s 2 4 > 


tz l 


vp>- p K» p ^> f^n 

[( s ,4> 


p 

V* 


] 


_ p S ? L 
- (S 2 +£) e 2 


A 2 (P) 


-k, ak 4 .v 2 ' 

’ plp+r-j Kp+r’ 2 ) 



(C-3) 
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where 


5,(p) - (C + R) , 

S 2 (p) = (C - R) , 



and r-j r 2 are given by Equation B-14. 
ii - Poles on the Complex Plane 

The straightforward process of finding the poles or "singular points" of 
the transformed temperatures T. (X,p) and T Q (X,p) and following it by application 
of the "residues" theorem, is the fastest procedure to solve for the inverse 
Laplace transformation. 

The values of p which result in an infinite value for Tj(X,p) or T Q (X,p) are 
herein called the poles and given by 

CD Pf=0, (2) p 2 =-r 1 , (3) p 3 =-r 2 , and C4) p 4 = p m 

where p m is an infinite set of complex poles that satisfy the equation 
Pm R m L p m -R L 

cc + r + R m> 6 s (C + v 1 - R m> e • 

V ""p^ pi ~ ~ . . . . (C-5 ) 

(ya+ocya + c + aic,) 

Since Tj ( X , p ) and T Q (X,p) are analytic everywhere in the region including the 
simple poles listed above, then the inverse transformation T(X,9) is given by the 
sum of residues of ^ p0 T(X,p)J at the poles of T(X,p). 
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Hi - Residues of the term [ A^(p) e + A 2 Cp) j : 


This term appears in Equations B-20 -21 for the image temperatures 

Tj and T 0 , respectively, It can be written as: 


r Si x SnXq 

L A 1 (p) e 1 + A 2 (p) e d J - 


2 f ( C+R ) X r p \\* 

AK 4 V^|e [(C+£ - R) e 


(C-R)L 


j] 


(C- 

e 




p(p+r 1 )(ptr 2 ){(C4^+R) ( e +R)L - (C+§- - K) < e" R)L } 


.(C-6) 


The cor resigning residue to each pole is listed as follows: 

(1 ) Residue at the Pole: p 3 0 

This represents the steady state value of this term and is given by 
- AK, ( C sinh R, (L-X) + R,cosh R, (L-X)l 

residup ss — -i ^ l w , 2 - 

(p=0) (C sinh R^L + R^osh R ] L) 

Where R-j = V C (C+2Ki ) 


. . . (C-7) 


(2) Residue at the pole : p = -r-j 

This is given by 

k,'ak 4 v 2 ; r i° 

residue = —7- — -V 
(p=-r-j) 

independent of the distance X. 


. . . (C-8) 


(3)» Residue at the pole : p = -r 2 

This is given by 

-k,ak 4 v 2 ; r 2 9 

residue = -—7- — ;n 

(p=-r 2 ) r 2^ r 2“ r l ' 

independent of the distance X. 


. . JX-9) 
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(4) Residues at the set of complex poles p m : 

This Is given by 

r V 1 ' 

p m (C+ P m + R ) ^ C+ *V L 
. m e 


S? 


p m 0 

r m 

K n AK, * e . R, 


m 


(C+R )X 
e 


in 


{l+ZLt^t C +1^)} (R 2 - C 2 ) 


which Is reduced to zero at the section X=d. 


1v - Residues of the term 






e + 




(c+fr - R) e 


This term appears In Equation B-21 for the image temperaturi 
written as: 


. m 4 v 


2 


(C-R)X 


(C+R)L 


e [2K 1 (C+^)e + f(C+f - R)1 - e ' UK^C+J) 


(C+R)X 


(C+R)L 


p(p+r-|)(p+r 2 ) • [ (C+y- + R) e - (C+y- - R) e 


(C-R)L 


The corresponding residue to each pole is listed as follows: 

(1 ) Residue at the pole p s 0 : 

This represents the steady state value of this term and 

cx 

AK a . e • si nh(L-X)R n 

residue = — 

(p=0) (C si nh R-jL + R^ cosh R-jL) 


where 

R 1 = Vc(C+ 2^) 

(2) Residue at the pole p - -r- | 
This is given by 

-^0 

-M 4 V e 1 

residue = 

(p= -r 1 ) ( r 2" r l^ 

independent of the distance X. 
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(C-R )X 
e m 

. . . (C-10) 

2*' 

a T 0 and is 

<C-R)L p# p v 
e + y{C+y +R) ] 

. . .(C-ll) 

is given by: 

. . . (C-12 ) 

. . . (C-13) 
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(3) Residue at the pole p « -r*; 


This Is given by 

-r 0 Q 

AK 4 V e 

residue = — 

(p= -r 2 ) (rg-r^ 

Independent of the distance X. 

(4) Residue at the set of complex poles p ^; 

This is given by 


v - Special Solution at Section X-Q 


. . . (C-14) 


p m 9 f p m \ c Jjn , n \~~i C+ V L I f (C+R )X p m 
^ AK 4 .R m e I (C V +, V e jje 

{R m 2 - C 2 )f 1 + 2L (^ + C + It,)} 


. . . (C-15) 


Combining Equations C-7, 8, 9, 10, 12, 13, 14 and 15 and substituting X c 0 


T-|(0,9) = T^ (0,0) = constant 


AK a tanh R,L 

T q (0,Q) = T (0,0) + —5 I + 


p 0 1 p 

m ? m 

AK 4 tanh R,l y R‘ (C-tfl+R^e 

(C tanh R,L + R ] ) m (R’jj - C 2 )jl +7L(^ia + C + K,JJ 

. . . (C-16) 


AK. V 

♦Inverse Laplace transform of )( 


Inverse Laplace transform of 


AK,V -r n 0 -r 9 0 
tr&) le • e > 


_ ak 4 AK 4 K,V 2 ; r l 9 ak 4 k,v 2 e ' r 2 0 
” 2 ) " 2C ■ r^rj-rp + r 2 (r 2 -r-J) 
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vi - Determination of the set of complex poles 

The infinite set of poles that satisfy equation C-5 can be found 
numerically by solving for the complex number Z such that 


si nh 2 Z = 


4L 2 K n 2 


Z = 2LR, 


m 


. . . (C-17) 


where Z = Z v + jZ is a complex number. The components Z and l are determined 
a y x y 

by solving numerically the following simultaneous equations: 


2LK-, sinh Z u cos Z t . - l = 0 
i x y X 

2LK-, cosh Z tf . sin Z - Z = 0 
i x y y 


. . . (C-18) 


The roots of Equation C-18 are only functions of the constants L and Kp and 
independent of the distance X, The sets p m and R^ can thus be followed from 
Equations (C-5) and (C-17). 
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Appendix D 

Inverse transformation from the frequency domain to the time domain for 


flow pattern (2) . 

i - Boundary and Initial Conditions 

With the image temperatures Tj (X,p) and T 0 (X,p) as determined from 
Appendix B, Equations B-19 and B-20 are subject to the following two conditions: 

(1 ) A boundary condition : 

At (X = L) and at all times 
T.|(L,e) = T q (L,Q) 
or 

T,a iP ) -t 0 o.,p) . . -cd-i> 

(2) An initial condition : 

At the inlet section (X=0) and at all times, the temperature T Q {0,9) is 
kept fixed, i .e. , 

T 0 (0>0) = T o (0,0) = constant 


or 


, T o (0,0) 

V».p) ■ -p — 


. . , CD-2) 

The arbitrary constants A-| Cp ) and A 2 Cp) for flow pattern. (2) are then given 
by: 


A 1 (p') = 


AK 4 V^ 


A 2 Cp ) = 


p(p+**1 )(p+r 2 l 


- MjjV 2 
pCp+r 1 )Cp+r- 2 > 


(C-*y -R)C y +K-| ) e - ^(Ki+C+^- -R) 

-CL r p Rl p W 
(C+y+R) e - (C+J- -R) e 


(C+y- +R) Cy- +K-| ) e 


CC+R)L 


^(K-j+C+y+R) 


where R = ~\j (y- +C } (^-+0 +2Kjj ) 


-CL 

d RL D -RL 1 


e 

(C+jJ+R) e - (C+£ -R) e 

J 


.(D-3) 
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11 - Poles on the complex plane : 

The values of p that lead to Infinite, values for Tj or T Q are identical 
to those found in Appendix C for flow pattern (1), 

I -(C+R)X -(C-R)X \ 

ill - Residues of the term l A-j(p) e + A 2 (p ) e I. 


This can be written as: 


9 -X(C+R) 
AK 4 VS e I 

n p (R-C)L p p 

[C-»y -R)( y +Ki ) e - y^K-j+C -ty -R) 


P(p+r 1 )(p+r 2 : 

p RL p -RL" 

) (C +£ +R) e - (C +y -R) e ; 

-CL 

e 


, (R-C)X 
AK 4 V . e 

p p -(C+R)L „ 

L (C+y +R)(f + K| ) e - y 

( K] + C + y + R) 

ptp+^Mp+rg) 

“ rC p 

(C+^ +R) e - {C+f -R) e 

-CL 

e 


(1 ) Residue at the pole p * 0 : 

This represents the steady state value of this term and is given by 

AK 4 -CX {c sinh R^L-X) - R 1 cosh R^L-X)} 

2C~ e (C sinh R-jL + R^ cosh R-|L) 

at X = 0 this is reduced to 


residue = 


AK 4 tanh R^L AK^ 


(C tanh R-jL + R-j J Zc" 

(2) Residue at the pole p = -r- j 

This is given by: 

AK 4 .K,V 1 2 -r-,0 
residue = — -r- — — -x e 
(p-r, ) r l lr 2 'V 

indejtehdent of the distance X. 


. . . (D-5) 


. . . (D-6) 
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